Universality of a family of Random Matrix Ensembles with logarithmic 

soft-confinement potentials 
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Recently we introduced a family of U (N) invariant Random Matrix Ensembles which is character- 
ized by a parameter A describing logarithmic soft-confinement potentials ^ [In H]^^^^\X > 0). 
We showed that we can study eigenvalue correlations of these "A-ensembles" based on the numer- 
ical construction of the corresponding orthogonal polynomials with respect to the weight function 
exp[— (In a;)^"'"'^]. In this work, we expand our previous work and show that: i) the eigenvalue density 
is given by a power-law of the form p{x) oc [Ina;]^"^/^; and ii) the two-level kernel has an anomalous 
structure, which is characteristic of the critical ensembles. We further show that the anomalous part, 
or the so-called "ghost-correlation peak", is controlled by the parameter A; decreasing A increases 
the anomaly. We also identify the two-level kernel of the A-ensembles in the semiclassical regime, 
which can be written in a sinh-kernel form with more general argument that reduces to that of the 
critical ensembles for A = 1. Finally, we discuss the universality of the A-ensembles, which includes 
Wigner-Dyson universality (A — >■ oo limit), the uncorrelated Poisson-like behavior (A — i> limit), 
and a critical behavior for all the intermediate A (0 < A < oo) in the semiclassical regime. We also 
comment on the implications of our results in the context of the localization-delocalization problems 
as well as the A'' dependence of the two-level kernel of the fat-tail random matrices. 

PACS numbers: 72.15.Rn, 05.40.-a, 05.60.-k, 05.90.-l-m, 71.30.1i 



I. INTRODUCTION 

Random matrix theory (RMT) deals with the statisti- 
cal properties of eigenvalues and eigenvectors of random 
matrices drawn from a certain probability measure. The 
theory was successfully applied by Wigner— in the 1950s 
to describe the spectral properties of complex many-body 
nuclei where the underlying Hamiltonians are so compli- 
cated that a useful way to study the system turned out to 
be through a statistical treatment of the Hamiltonians^. 
Since then, RMT has been applied to a wide variety 
of systems in diverse areas, including e.g. many-body 
atoms and nuclei, quantum chaos, mesoscopic disordered 
conductors, 2-D quantum gravity, conformal field theory, 
chiral phase transitions as well as zeros of Riemann zeta 
function, scale-free networks, biological networks, com- 
munication systems and financial markets^i^. This broad 
range of applicability of RMT in seemingly unrelated ar- 
eas highlights the universal features of the correlations 
of the eigenvalues in RMT. Within the classical Gaus- 
sian model pioneered by Wigner, these correlations are 
known as the Wigner-Dyson (WD) statistics of the Gaus- 
sian ensembles, which are qualitatively different from the 
statistical features of completely uncorrelated eigenvalues 
given by the Poisson statistics. 

In the past decades, many attempts have been made 
to construct generalized random matrix ensembles that 
incorporate power-law or fat-tail distributionsi^"— . The 
significance of such generalization beyond the Gaussian 
ensembles is mainly two fold. First of all, it is conceiv- 
able that the wide applicability of classical or Gaussian 
RM models is closely linked to the prevalence of nor- 
mal or Gaussian distributions in nature, a consequence 
of the central limit theorem^. It has, however, not been 
fully investigated if there is a counterpart of the Gaussian 



ensemble as implied by the structure of the generalized 
central limit theorem^, comparing e.g, the Gaussian and 
the Levy basins. Second, there have been expanding in- 
terests in RMT applications to generic complex systems 
such as financial markets, scale- free network, earthquakes 
etci^i^ that feature fat-tail noise, but the relevance of the 
classical RM models in these systems seems questionable. 
This is because the Gaussian ensemble is based on the as- 
sumption that the systems are characterized by Gaussian 
noise, which is clearly not suitable for systems with fat- 
tail noise where the occurrence of extreme events are not 
as rare as expected from normal distributions. So far, fat- 
tail distributions in random matrix ensembles have been 
carefully incorporated in some limited cases and the cal- 
culation of the correlation functions of eigenvalues have 
been carried out for certain special cases!. However, the 
question regarding the universality of the correlations of 
the eigenvalues remains unresolved. For Gaussian ensem- 
bles it is the well known two-level sine kernel that estab- 
lishes the universality of the correlations in the properly 
scaled large N (matrix size) limit; it is not clear if there 
exists a similar universal two-level kernel for the power- 
law or the fat-tail ensembles as well. 

In fact, the suggestion of a novel universality beyond 
the Gaussian ensembles comes from the study of the An- 
derson transition in the disordered electronic systems^. 
In these systems, the Gaussian ensemble is only relevant 
in the metallic regime where all the eigenstates are ex- 
tended across the entire system and correlations of the 
corresponding eigenvalues are well described by the WD 
statistics. As the disorder is made strong enough, the 
eigenstates become localized and thus the eigenvalues 
become uncorrelated. Especially at the delocalization- 
localization transition, it has been established that the 
correlations of the eigenvectors exhibit novel features^ 
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such as multi-fractality and the correlations of the eigen- 
values lead to a level compressibility that is intermedi- 
ate between WD and Poisson statistics. Similarly in the 
studies of quantum chaos, energy level statistics of sys- 
tems that are intermediate between chaotic and regular 
states also require generalization beyond WD and Pois- 
son statistics^. In these contexts, extensive studies have 
been carried out to construct a parametric generalization 
of RM models that cross over from WD to Poisson^ii^ as a 
function of the parameter. Some of these generalizations 
indeed capture the essential features of the critical statis- 
tics, among which the family of q-RMEsiiii^ provides a 
particularly valuable insight. Within the common frame- 
work of rotationally invariant RM models^^ the q-RMEs 
show how the universality of the Gaussian ensemble char- 
acterized by the well-known sine kernel breaks down and 
eventually gives rise to a different kind of universality for 
the critical ensembles, characterized by a one-parameter 
two-level sink kernel. In particular, the rotationally in- 
variant RM models are characterized by a "confining po- 
tential" which defines the weight function of a set of 
orthogonal polynomials; the key difference between the 
Gaussian and the critical ensembles comes from the fact 
that the corresponding orthogonal polynomials, namely 
classical vs. g-orthogonal polynomials, respectively, pos- 
sess qualitatively different asymptotic propertiesi^ii^. 

In the motivation to investigate the universality as- 
sociated with fat-tail or power-law RMEs, we intro- 
duced a family of U {N) invariant random matrix ensem- 
bles characterized by an asymptotic logarithmic poten- 
tial V{H) = A\[nxY+^ with A > Oi^, which we wih re- 
fer to as "A-ensembles"— . The reason for such a choice 
of the potential is based on the following few observa- 
tions. First, it is known that for V{H) cx [Ini?]^ corre- 
sponding to A = 1 (the critical ensemble), the eigenvalue 
spectrum is given by inverse power-law distribution, as 
known by the mean-field theoretic approacbii. Second, 
for V{H) c:^ NlnH corresponding to the A — >■ limit with 
the constant A being order of N (free Levy matrices), 
the spectral density is given by the fat-tail distributions. 
Third, in the limit A 3> 1, it is expected that the con- 
finement potential may grow sufficiently strong, thereby 
approaching the Gaussian limit. The fact that such a 
parametric generalization connects various existing RM 
models is interesting since the model allows us to explore 
any possible novel universality associated with fat-tail 
RMEs and with logarithmic soft-confinement potentials 
within the rotationally invariant RMT framework. 

It is well known from the mean-field theoretic 
approacbii that RMEs with logarithmic soft-confinement 
arc characteristically different from those with strong 
confinement potentials given by V{x) = a > 1. The 
reason is that for the soft-confinment potential, the eigen- 
value density does not depend on N and has a non-trivial 
functional form, which is not translationally invariant. 
It means that the required unfolding is non-trivial. In 
fact, the non-trivial unfolding procedure give rise to the 
deviation from the WD statistics. In the following, we 



will show by studying the two-level kernel that the A- 
ensembles exhibits interesting deviations from WD statis- 
tics. Especially, the presence of the anomalous compo- 
nent in the two-level kernel is a common characteristics 
of all the logarithmic confinement potentials. 

The paper is organized as follows. In sec II wc provide 
a brief review of the orthogonal polynomial method and 
in sec III wc discuss the model ensemble. In sec IV we 
will show the results and finally in sec V wc will discuss 
the results with concluding remarks. 

II. ORTHOGONAL POLYNOMIAL METHOD 

We consider the set of U{N) invariant Hermitian ma- 
trices H with the following probability measure 

PN{H)dH (xe-^'^^''"^^dH, (1) 

where the confining potential V(x) is a suitably increas- 
ing function of x, tr is the matrix trace and dH the Haar 
measure. In the eigenvector basis, the joint probability 
distribution of the eigenvalues {xi}{i = 1,2,... N) of the 
matrices can be written in the formic 

N N 

PN{{x^})ocY[{x,~x,fY[e-''^-'\ (2) 

i<j 4=1 

Here the factor Yii'-'^i ~ ^j) ^^'^ Vandermondc de- 
terminant. Given a set of (monic) polynomials Pn{x) 
that are orthogonal with respect to the weight function 
wi^x) = e~^(^^, i.e. 

oo 

j e~^^'''>Pn{x)pm{x)dx = Smn, (3) 
— oo 

the n- level correlation of eigenvalues Tln{xi,X2, ■■■■,Xn) 
can be written in a compact form in terms of the the 
two- level kernel in the following manner, 

Tlnixi,X2,—;Xn) = dct[K N {xi , Xj)]{i,j=o „}. (4) 

Here the two-level kernel Kn{x, y) is defined as 

KN(x,y) = ^ 1pn{x)llJn{y), (5) 
n=0 

where the "wave function" ipn{x) = Pn{x)e^^^^^^'^ . By 
using the ChristofTel-Darboux formula, the kernel is sim- 
plified to (for the monic polynomials) 

, . ipN{x)ipN-i{y) - tpNiy^N-iix) 

KN[x,y) = . (6) 

x-y 

In general, the large N asymptotic behavior of the or- 
thogonal polynomials with respect to any weight func- 
tion characterized by V{x) ~ x", a > 1, have behavior 
qualitatively similar to the Hermite polynomials^, such 
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that the asymptotic behavior of the wave function in the 
iV — >■ oo hmit is given by 



Tp2N{u) ^ cos(7ru) ; tp2N-iiv) sin(7rw) 



(7) 



where u and v are scaUng variables, namely, du = 
Kao{x,x)dx and dv = Kooiy,y)dy such that the mean 
density Koo{u,u) is unity in the N — >■ oo limit. In this 
limit, the two-level kernel in the scaled variables become 



sin{7rs) 



TTS 



{u — v) 



(8) 



which is the celebrated sine kernel. Thus, in the frame- 
work of orthogonal polynomial method, the universality 
of random matrix ensembles is traced to the similarity 
of the scaling behavior in the large limit of the poly- 
nomials corresponding to the probability measure under 
considerationi^. In particular, the orthogonal polyno- 
mials corresponding to all 'Freud-type' weight functions 
e~^^^\ with monotonically increasing polynomial V{x), 
share similar asymptotic behavior in the large TV limit. 
Thus, they give rise to the same correlations as the Gaus- 
sian ensembles. 



III. THE A-ENSEMBLES 

The generic choice of the confining potential V{x) 
that gives asymptotic logarithmic behavior is V{x) = 
A[\nx]^'^'^ . However, it has an unphysical singularity 
at the origin so that we need to regularize it in a cer- 
tain way. One possible way to do it is choosing, e.g., 
V{x) = A[ln(l -I- 3::)]^+'*' but there are a variety of other 
forms that differ by the regularization behavior in the 
vicinity of origin, which will not change the characteris- 
tics of the A-ensembles. For our study, we particularly 
choose the following form of potentiaU^: 

V{x) = i[sinh"^a;]i+^; A > 0, 7>0 (9) 

7 

The merit of choosing Eq. (jH]) is that for A = 1, it co- 
incides with the one possible form of the weight func- 
tion of the g-RMEsiSS so that we can compare our re- 
sults with those of (/-RMEs. For the g-RMEs, the mathe- 
matical properties of the corresponding orthogonal poly- 
nomials, "the Ismail-Masson g-polynomials"— are well 
known, which leads to the two-level sinh kernel in the 
limit 7 ^ 27r, given by 



7 sin(7rs) 
27rsinh(7s/2)' 



(10) 



As is well known, this recovers the Gaussian sine kernel 
in the limit 7 = 0. 

For A 7^ 1, the corresponding orthogonal polyno- 
mials, which we refer to as "A-generalization of q- 
polynomials"— are known via the recurrence relation of 
the orthogonal polynomial, i.e., 

x(l)nix) = 0„+i(a;) -I- S'„0„(a;) -f- Rn(f>n~i{x), (11) 



where Sn and i?„ are real. For symmetric weight func- 
tions all Sn are zero, and the polynomials are determined 
by the coefhcient The following shows comparison 
among the Hermite polynomials, g-polynomials, and the 
A-generalization of (/-polynomials: 



Rn oc n Hermite polynomials. 



oc e 



7n 



(12) 
(13) 

A-gcneralization of (7-polynomials.(14) 



cx e'''" q-polynomials, 



We can see clear distinction between i?„ of the Hermite 
polynomials (linear in n) and that of the A-generalization 
of g-polynomials (exponential in n^/^) which reduces to 
that of g-polynomials for A = 1. The significance of i?„ 
is that it determines the upper bound of the spectral 
density and thus the scaling behavior of the bulk of the 
spectrum in the large N limiti^ii^. For example, for the 
Gaussian ensembles characterized by V{x) = x^ (Her- 
mite polynomials), the upper bound of the spectral den- 
sity is ^/N in the large N limit. Thus the normalization 
condition of the spectral density requires the bulk of the 
spectrum to grow as \/N as in the semi-circle law. On 
the other hand for the logarithmic soft-confinement po- 
tentials, the spectral edge grows at an exponential rate 
e^"^^^ . The bulk of the spectrum does not scale as N . For 
example, for V{x) oc [Inx]^, the spectral edge grows as 
V and the bulk of the spectrum is given by p{x) (x 1/x 
which does not depend on N . 



IV. RESULTS 

The calculations of the density and the cluster function 
arc performed based on the numerical constructioufi^ of 
the A-gencralization of the g-polynomials. For instance, 
the density p{x) = K]\j{x, x) can be obtained by summa- 
tions of products of the wave function as in Eq. ([5]). The 
two-level cluster function 



Y{u,v)=[Kiu,v)f 



(15) 



can be obtained based on Eq. ([5|) with the unfolding map 
u{x) = dxp{x). 



A. Eigenvalue Density 

In Ref. [l^ it was observed that the density of the A- 
ensembles is given by a pure power-law, e.g., p{x) = . 
For A = 1, e = and for A > 1 and A < 1, > o'^and 
6 < 0, respectively. However, careful further investiga- 
tion shows that this observation is only approximate and 
a more accurate form of the eigenvalue density is given 
by 



p{x) oc for X > A 



(16) 
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FIG. 1: Eigenvalue density for A < 1, namely A = 0.9 (7 - 
0.75), 0.7 (7 = 0.50), 0.5 (7 = 0.25) as well as 1 (7 = 0.75). 



FIG. 3: Logarithmic dependence of the eigenvalue density for 
A-ensembles for A < 1. 
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FIG. 2: Eigenvalue density for A > 1, namely A = 1.1 (7 = 
0.75), 1.3 (7 = 2.00), 1.5 (7 = 4.00) as well as 1 (7 = 0.75). 



FIG. 4: Logarithmic dependence of the eigenvalue density for 
A-ensembles for A > 1. 



where the lower cutoff A depends on the parameters in 
y{x)- For example for A = 1, A = 7 and for A ^ 1, A 
is some function of both 7 and A, which depends on the 
specific regularization of Yix) at the origin. 

The above form is suggested by the mean-field 
approacbi^i for A > 1. The validity of this form for all 
A > can be checked by considering the normalization 
condition of the spectral density 



2 / p[x)dx = N. 
'0 



(17) 



Here the factor 2 comes from the fact that p{x) is sym- 
metric around origin. The upper bound Dn is given by 
the largest zero of the orthogonal polynomials of order N, 
namely Dn oc \/Rn- As pointed out, i?„ oc exp[n-^/'^]. 
We notice that Rcf. [23| studied the largest zeros of 
the orthogonal polynomials to the weight function of 
exp[— c(lna:)™] for c > and m a positive even integer 

and showed that the largest zero is of order exp(n™-i ), 
which is the same behavior as the coefficient i?„ of the 
A-generalization of g-polyonomials in Eq. (|14p . Thus, our 



results seem to imply that the results of Ref. [20[ can be 
extended to an arbitrary real A > 0. 

In order to test the validity of Eq. numerically, we 
assume that the density is in general given by 



p{x) 



X + A 



where /(x; A) is a slowly varying function of x. 
xp{x) will have the form 



xp{x) = 



X 



x + A 



/(•t;A), 

X > A. 



(18) 
Then 

(19) 
(20) 



Thus, the function /(x; A) is identifiable in the large x ^ 
A regime oixp{x) vs. x plot. Figure [T] shows the logarith- 
mic behavior of xp{x) for A < 1 (A = 0.5, 0.7, 0.9, and 1), 
while Figure[2]show it for A > 1 (A = 1.1, 1.3, 1.5, and 1), 
respectively. For all the cases, we chose 7 ~ 0(1), which 
ensures the cut-off A = 0{1). 

To further investigate if f{x;X) oc [Inx]'*^"^ for large 
X {x ^ 1), we plot \n[xp{x)] vs. Inlnx and fit it in the 
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FIG. 5: Anomalous component of the cluster function for A- 
ensembles for A < 1 (7 = 0.50). 




FIG. 6: Anomalous component of the cluster function for A- 
ensembles for A > 1 (7 = 1.50). 



range 10 < x < lO"'. Figures [3] and |4] show the expected 
linear behavior, establishing the validity of Eq. ()16[) . 



B. Ghost correlation peak 

For A =1, the existence of a "ghost correlation peak" is 
well knowni^iiiii^iiS^. The ghost peak was first discussed 
in Ref. 22| though it was contained implicitly in the 
exact result of Ref. [T^]. While the normal part of the 
two-level cluster function for the critical ensembles for 
uv > and |w — w| << it is given, in the 7 27r^ limit, 
by the sinh kernel 



7 sm[7r(u - 



27r sinh[^(M — v)] 



(21) 



there exists an anomalous part of the cluster function for 
uv < given by 



y:{u,v) 



' 7 sin[7r(M — v)] 



27r cosh[^(M + v)] 



(22) 



In fact, the presence of such long range correlation is 
required by the normalization sum rule^^ 

/oo 
du[Y:i{u,u')+Y,^u,u% (23) 
00 

The deficiency of the sum rule 

/oo poo 
duY,1{u,u')^ / duY^{u,u') (24) 
-00 J —00 

is related to certain characteristics of the critical 
statistics^: i) the level compressibility in the number vari- 
ance S(L) within a range L and ii) the mult i-fr act ality 
of eigenvectors. In particular. 



X 



d^{{L)) ^ d-D2 
d{L) 



2d 



(25) 



where d is the system dimension and the fractal dimen- 
sionality Dp determines the scaling behavior of the mo- 
ments of the inverse participation ratio via 



(26) 



It turns out that A-ensembles possess such nor- 
mal/anomalous structure for A 7^ 1 as well. Figures[5]and 
inishow the numerical evaluation of 1- Y^{u, v) {u > and 
V < 0) for a symmetric range around v — for varying 
A values for a fixed 7 (7 = 0.5 in Fig. [S] and 7 = 1.5 in 
Fig. [6]). As the figures show clearly, the magnitude of the 
ghost peak depends on A in a significant way; for A < 1 
(Fig. [5]), the peak is more pronounced than that for A = 1 
and for A > 1 (Fig. ^ it is the opposite. 

The observation that such long range correlation lead- 
ing to the ghost peak is preserved for all A 7^ 1 suggests 
that such features are common to all logarithmic con- 
finement potentials. In other words, once the critical 
ensembles break the U (N) symmetry^ of the Gaussian 
ensembles with the introduction of the parameter 7 (or 
q = e^'''), the A ensembles remain in this broken symme- 
try family. The fact that as A becomes large, the ghost 
peak shrinks seems to imply that the U {N) symmetry 
might become fully restored in the limit of A — cxj. This 
expectation seems consistent with the asymptotic behav- 
ior of two-level correlation in the limit of A — > cx) that will 
be shown later in section D. 

The A-ensembles are all "critical" in the sense that 
the normal part of the two-level kernel violates the sum 
rule that can be associated with the characteristics of 
the critical statistics such as level compressibility and 
multi-fractality. In particular, the fact that the violation 
of the sum rule is controlled by the parameter A, e.g.. 
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Gaussian 
?L=0.9 




FIG. 7: Normal component of the cluster function for A- 
ensembles for A < 1 (7 = 0.50) 




FIG. 8: Normal component of the cluster function for A- 
ensembles for A > 1 (7 = 1.50) 



< xW < 1 is intriguing. As mentioned above, A seems 
related to the degree of the U{N) symmetry breaking 
and thus, indicative of the non-trivial character of the 
eigenvector correlations, namely the multi-fractal dimen- 
sionality as can be seen immediately from Eq. (|25p and 
the A-dependent sum-rule deficiency. In this regard, the 
study of the dimensional dependence of the critical statis- 
tics will be important to further understand the role of 
the parameter A since the multi-fractal dimension of the 
eigenvector correlations at the critical states is dependent 
on the dimensionality of the system^^. 

Figures [7] and [8] show the comparison between the nor- 
mal part of cluster function for A 7^ 1 and that for A = 1 
as well as that of the Gaussian ensemble. We note that 
the nodes of the cluster functions occur at integer val- 
ues on the X axis for all A values as well as the Gaussian 
case, which is the indication of the semi-classical behav- 
ior. The peak height and position between the nodes 
gradually decrease and shifts as A decrease for a given 7 



value (7 = 0.5 in Fig. [71 an d 7 
already shown in Ref. [iq . 



1.5 in Fig. [8]), which is 



C. Identification of two- level kernel 

As mentioned before, for the logarithmic soft- 
confinment potential, the eigenvalue density is a non- 
trival function that do not depend on N as shown in 
section A. The deformation from the Gaussian universal- 
ity, the 'sine kernel', is expected due to the non-trivial 
unfolding procedur o^^i^^ . In the following, we will iden- 
tify the two-level kernel of the A-ensembles in the semi- 
classical regime and show that the proposed form is in 
good agreement with the numerically-obtained cluster 
function (thus the two- level kernel). 

In the semi-classical regime (7 < 2tt for the critical 
ensembles corresponding to A = 1), the kernel for an 
arbitrary weight function can be written a o^^'^^ 



K{u,v) 



sin(7r(M — v)) 
x{u) - y{v) 



(27) 



For example, for the critical ensembles (A = 1) the un- 
folding variable, in which the density becomes uniform 
and unity, is given by 



Mo 



t 



-dt = lna;/a;o. 



(28) 



This leads to x — x^e^ and hence to the sinh kernel. 
In a similar manner, we define the unfolding variable for 
the A-ensembles as 



cflnt 



A-l 



t 



-dt 



A 



([ln.T]^-[lnxo]^l29) 



where uq — p[x)dx and Xq 3> A, the cut off of the 
density near the origin. The constant c cx p{x ~ 0). 
Note that p{x = 0) is a hmction of the parameters that 
do not depend on N . For the critical ensembles, p(0) — ^ 

and for the A-cnscmbles, p(0) = g{X^^) where g is some 
function of A and 7. Rewriting the original variable x in 
terms of the unfolding variable, we obtain 



exp "f*{u' + af/^ 



(30) 



where 7* 



,1/A 



and the a = ■|[lna;o]'''- For A = 1, it 



reduces to the unfolding for critical ensembles. Eq. (P7)) 
then suggests the two-level kernel for A-ensembles to be 
of the form 



r sin[7r(ft — v)] 
27r sinh[E-({t — v)] 



(31) 



where we have defined ii ~ u' -\-a and v = u' + a. and the 
parameter F is introduced to satisfy the condition that 
(t{u) = K^{u,u) = 1. It then follows that 



T{u,v) = 7* 



,71/A ,71/A 



(32) 
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Thus we arrive at the regular component of the two-level 
kernel of the A-ensembles 



r(u, u) sin[7r(u — u)] 
27r smh[^^{il-v)y 



uv > 0, (33) 



while the anomalous component of the two-level kernel 
of the A-enscmbles is given by 



T{u,v) sin[7r({t — w)] 



cosh[^^iu + v 



— ; uv < 0. (34) 



We note that the kernel Eq. ((33|) reduces to the well- 
known sinh kernel for A = 1, where r(u, 5) — >■ 7* = 
7. In general it is not translationally invariant, e.g. for 
A = 0.5,r(w,v) = -f*{u + v) and for A = 2.0,r(u,S) = 
—J- — 7^. However, if we choose -D to be a fixed value, i.e. 

Vu-\-Vv 

v = v^uo + a = a, which is the same as choosing u = uq, 
then u = u — UQ + a = u — v + a = s + a. In this way, 
the function r(u, v) and g(u, v) can be written in terms 
of a difference variable s = u — v = u — v alone with a 
constant v = a that serves as a fixed reference point. 



r(s + a,a) — r(s, a) = 7' 



(s + a) 



1A_„1/A 



Then 



r(s,a) sin(7rs) 



2tt 



sinh[ 



(35) 



(36) 



Figures l9l and fTOl show the fitting results with the clus- 
ter function corresponding to the kernel Eq. 1311 for dif- 
ferent A values. They show a fit with the numerically- 
obtained two-level kernel, from the A-generalization of 
the g-polynomials and Eq. ([5|), with fit values 7* ~ 7 = 
0.5 and a ~ 2 for A < 1 (Fig. O and 7* ~ 7 = 1.5 and 
a ~ 2 for A > 1 (Fig. In Figures [IT] and [111 we show 
that in similar range of the parameters, the anomalous 
components of the kernel also agree very well with the 
proposed form of the kernel Eq 



D. Novel universality 

In the following, we will discuss novel asymptotic be- 
havior of the two-level kernel of the A-ensembles by fur- 
ther examining Eq. ([36]). First, we note that 



T[s,a)(xs>- ^ for s 3> a 

leads to 

Kn{s, a) (X s^^^e^'^ ^ for s S> a. 
Looking at the A-dependcnce, we observe that 



(37) 



(38) 



K„{s,a) (X e 
cx e~ 



for A 
for A = 1 



Gaussian 
?.=1.0 
X=0.9 (fit) 
X=0.8 tit) 
1^0.7 (tit) 




FIG. 9: Fitting results for normal component of the cluster 
function for A-ensembles for A < 1 (7 = 0.50). The lines 
are obtained from Eq. (jSj using the A-generalization of the q- 
polynomials, while the points are obtained from the analytic 
form of the cluster function given by Eq. (|33[l. 




1 1.5 



FIG. 10: Fitting results for normal component of the cluster 
function for A-ensembles for A > 1 (7 = 1.50). The lines 
are obtained from Eq. ^ using the A-generalization of the q- 
polynomials, while the points are obtained from the analytic 
form of the cluster function given by Eq. (|33[) . 



for A — > 00. 



(39) 



Thus, in the limit A — >■ 00, we get back the well-known 
decay of the Gaussian ensembles (WD universality). This 
is consistent with the fact that the magnitude of the 
ghost correlation peaks become smaller as A increases, 
presumably disappearing in the Gaussian limit of very 
large A. For A — )■ limit, the asymptotic tail is given 
by the infinitely fast exponential decay, which suppresses 
the eigenvalue correlations, thereby leading to uncorre- 
lated Poisson-like behavior (Poisson statistics). In be- 
tween these limits, the large s behavior is governed by 
^ , which is a novel feature of A-ensemblcs. 
In a similar fashion, we can study the two-level density- 
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FIG. 11: Fitting results for anomalous component of the clus- 
ter function for A-ensembles for A < 1 (7 = 0.50). The lines 
are obtained from Eq. ^ using the A-generalization of the q- 
polynomials, while the points are obtained from the analytic 
form of the cluster function of Eq. (f34|) . 




FIG. 12: Fitting results for anomalous component of the clus- 
ter function for A-ensembles for A > 1 (7 = 1.50). The lines 
are obtained from Eq. ([5]) using the A-generalization of the q- 
polynomials, while the points are obtained from the analytic 
form of the cluster function of Eq. (|34p . 



(40) 



density correlation function 

R{u, v) = 6(u — v) — Y{u, v) 



In terms of the difference variable s = u — v, the two- level 
correlation function for the A-ensembles can be written 
as 



(41) 



For s ^ a, i?'^(s) ~ for all values of A, which is an 
expected feature due to the unitary symmetry of the en- 
sembles. For large s ^ a, the results are simply obtained 
from Eq. ^ 



FIG. 13: Number variance for A-ensembles for A =0.5, 0.8, 
1.0, 1.5, 2.0 (7 = 0.50). 



cx e for A 
cx for A 



(42) 



We can see that the asymptotic behavior of the two-level 
correlation function interpolates that of Gaussian ensem- 
bles for A — 00 limit and uncorrelated Poisson-like be- 
havior for A — > limit. This suggests that A-ensembles 
belong to the critical ensembles characterized by A which 
determines the rate of asymptotic exponential decay of 
the two-level kernel. 



E. Number variance 

As an example of the spectral measures obtained from 
the two-level kernel, the evaluation of the number vari- 
ance S(i) within a range L is shown in Figure [T51 It 
clearly shows that for all values of A, the number vari- 
ance is linear in L for large L. As A becomes smaller, the 
number variance shifts towards the uncorrelated Poisson 
distribution. This is consistent with the fact that defi- 
ciency of sum rule rule x(A) increases as A decreases as 
shown in section B, since the slope of number variance is 
directly related to the deficiency of the sum rule: 



X(A) 



d{L) ■ 



(43) 



R (s,a) cx e 



for A ^ 



V. SUMMARY AND DISCUSSION 

In this work, we further study an invariant-class of 
random matrix ensembles characterized by the asymp- 
totic logarithmic soft- confinement potentials introduced 
in Ref. [l^ which we refer to as the A-ensembles. As 
a first step, we carefully reinvestigate the spectral den- 
sity of the A-ensembles and show that instead of a A- 
dependcnt power law in a restricted regime, the spec- 
tral density is given more accurately by a power-law of 
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the form p{x) oc [lna;]'^~^/x. This result is suggested by 
the mean-filed approach and it can be checked by the 
normalization condition of the spectral density, the edge 
of which is determined by the coefficient i?„ of the A- 
generalization of the g-polynomials. 

Second, we show that the two-level kernel of the A- 
ensembles has normal/anomalous structure, which is 
characteristic of the critical ensembles. The anomalous 
component arising due to the sum rule violation is de- 
pendent on the parameter A for a given value of 7 which 
is kept 0(1) in the numerics; As the value of the A de- 
creases, the deficiency of the sum rule becomes larger. 

Third, we identify the normal and anomalous compo- 
nents of the two-level kernel in the semi-classical regime, 
which are given by Eq. ([33| and Eq. (p4|) that reduce 
to those of the critical ensembles for A = 1. Further, 
we show that the two-level kernel of the A-cnscmbles ex- 
hibits a distinct universal asymptotic behavior, shown 
by Eqs. ([5^ . which includes the Gaussian ensembles 
(A — ?> 00 limit), the critical ensembles (A = 1) as well 
as the uncorrelated Poisson-like behavior (A —J- limit). 
In particular, the large s behavior of the two-level ker- 
nel is governed by exp[— s^^"^], which is the characteris- 
tic feature of the A-ensembles. It is expected that the 
asymptotic tail of the spacing distribution is also given 
by a similar exponential form. 

Lastly, we show that the number variance is linear in 
L for large L. The slope of number variance or the level 
compressibility is dependent on the parameter A. As A 
decreases, the slope increases, which is consistent with 
the fact that anomaly or the deficiency of the sum rule 
increases as A decreases. 

All the features of the A-ensembles shown above such 
as the ghost peak, the deficiency of the sum-rule, the 
finite compressibility as well as the asymptotic expo- 
nential decay of the kernel seem to suggest that the 
A-ensembles belong to the critical ensembles character- 
ized by A. Therefore, it is expected that the A-ensemblcs 
would be relevant to the description of the critical states 
of the localization-delocalization problems in disordered 
systems. Since the critical level statistics are universal. 



depending only on the critical exponent and the dimen- 
sionality of system for a given symmetry class, it is con- 
ceivable that the parameter A can be associated with 
these parameters. In particular, the study of the dimen- 
sional dependence of the critical statistics is interesting 
since the parameters of the critical level statistics such 
as the critical exponent and the multi-fractal dimension 
are dependent on the spatial dimension of the system. 

At the same time, it is also interesting to see if A- 
ensembles are also applicable in the study of quantum 
chaos. Recently the critical statistics has been found rel- 
evant in some cases of quantum chaos as well^^. It turns 
out that the two level kernel of chaotic systems with log- 
arithmic singularity^ have the exact same form as that 
of the critical ensembles. 

Another important implication of the A-cnscmbles is 
in regard to the question that we posed earlier, namely 
if there is a universality of the correlations of the eigen- 
values associated with fat-tail distributions. Our results 
seem to suggest that a non-trivial A^-dependence of two- 
level kernel of the fat-tail RMEs is a generic feature 
within the framework of rotationally-invariant RMT. In 
the A — > limit of the A-ensembles, 7 is required to be an 
A^-dependent parameter to have the probability measure 
to be normalizable, which is the case for the free Levy 
matrices. Thus, the A^-dcpendence of the two-level ker- 
nel in this limit can be understood as a consequence of 
the presence of the A^-dependent parameter in the model 
ensemble, which cannot be simply scaled out. 

In conclusion, we show that the family of RMEs with 
the asymptotic logarithmic soft- confinement potentials 
characterized by A, called the "A-ensembles", connects 
the WD universality (A 00), the uncorrelated Poisson- 
like behavior (A — >■ 0) and exhibit a critical behavior for 
all the intermediate A value (0 < A < 00) in the semi- 
classical regime. We expect that further study of the A- 
ensembles and the A-generalization of g-polynomials will 
lead to a deeper understanding of the universality of ran- 
dom matrix ensembles in general. 

We gratefilly acknowledge valuable discussions with 
V.E. Kravtsov. 
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